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ABSTRACT 

Computational results are presented for the 
performance and flow behavior of various injector 
geometries employed in transverse injection into a non- 
reacting Mach 1.2 flow. 3-D Reynolds-Averaged 
Navier Stokes (RANS) results are obtained for the 
various injector geometries using the Wind code with 
the Mentor’s Shear Stress Transport turbulence model 
in both single and multi-species modes. Computed 
results for the injector mixing, penetration, and induced 
wall forces are presented. In the case of rectangular 
injectors, those longer in the direction of the freestream 
flow are predicted to generate the most mixing and 
penetration of the injector flow into the primary stream. 
These injectors are also predicted to provide the largest 
discharge coefficients and induced wall forces. Minor 
performance differences are indicated among diamond, 
circle, and square orifices. Grid sensitivity study results 
are presented which indicate consistent qualitative 
trends in the injector performance comparisons with 
increasing grid fineness. 

INTRODUCTION 

Motivation 

Fluid injection into a supersonic crossflow has recently 
been a focus of investigations for application to 
ramjet/scramjet combustion. Quick penetration and 
mixing of injection fluid are metrics of interest here due 
to the short combustor residence times involved. 

There is also significant interest in the applications, 
albeit with lower speed supersonic crossflows (M < 2), 
to jet engine noise damping and thrust-vectoring. The 
interest in noise damping stems from the success of 
chevron nozzles, which have been shown to 
significantly reduce jet engine noise levels. Apparently, 
the large-scale streamwise vortices created by the 
chevroned-nozzle geometry induce jet exhaust- 
freestream mixing, which promotes noise reduction. It 
is thought that creating strong streamwise vortical 
structures (and the resulting increased mixing) via 
transverse injection near the trailing edge of a nozzle 
could promote noise reduction in an analogous manner. 


The application to thrust-vectoring within a nozzle by 
using an injected flow to induce a wall force is also of 
interest relative to improving nozzle performance or 
vehicle stability and control. 

Clearly, it is important to obtain a qualitative 
understanding of how the aforementioned performance 
characteristics (i.e., penetration, mixing, and induced 
wall force) are affected by injector configuration 
parameters. There are a wide variety of possible 
injector configurations including variations in orifice 
geometry, orientation relative to the primary flow, 
injector-to-freestream pressure ratios, etc. 

Literature Survey 

Tomioka, et al. [1], have experimentally investigated 
the sensitivity of penetration and mixing to injector 
configuration for a Mach 3 crossflow. The parametric 
study observed the effects of variations in jet-to- 
freestream dynamic pressure ratio, injector orifice 
shape, and injection angle. They noted better 
penetration for diamond-shaped orifice injectors versus 
equivalent area, circular orifice injectors. At large jet- 
to-freestream dynamic pressure ratios, the plume 
penetration of the diamond-shaped injector is 
comparable to that of the circular injector. 

Santiago and Dutton [2] experimentally studied the 
streamwise counter-rotating crossflow vortices 

generated during transverse injection into a supersonic 
crossflow. They assert that these vortices are the 
primary source of entrainment of crossflow mass into 
the injected fluid jet plume. Ben-Yakar et al. [3] found 
similar flow phenomena (e.g., counter-rotating vortices) 
for subsonic injection. 

Research Objectives 

The present study is focused on understanding the 
effects of variations of orifice geometry on the 
performance of flush-mounted injectors for a low 
supersonic flow (Mach = 1.2). The performance criteria 
were penetration, mixing, and induced wall forces. The 
injector-to-freestream static pressure ratio is roughly 
2:1. A numerical methodology is also developed to 
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provide a reasonably efficient and fair comparison of 
injector performance among the various candidate 
orifice shapes. Note that experimental validation of the 
numerical findings is not available at the present time. 

METHODOLOGY 
Model Geometry 

The model geometry (Fig. 1) consists of a Mach 1.2 
primary stream duct flow, and a subsonic (-Mach 0.3) 
secondary injector flow. The injector domains are 
connected to the duct bottom wall, and have uniform 
cross-sectional area (as shown). The duct domain is an 
8 by 8 in. cross-section, which extends for 96 in. All 
orifice geometry areas are 8 in. 2 . The choice of the 
relative size of the injector areas to duct width is 
arbitrary. Note that the sidewalls model neighboring 
rows of injectors (i.e., mirror images). The injector 
orifice geometries, depicted in Fig. 2, include circular, 
square, diamond, and eight slots with different length- 
to-width (L/W) aspect ratios. For example, the 1:8 slot 
(or latitudinal slot) is 1 by 8 in., and oriented across the 
entire width of duct, while the 8:1 slot (or longitudinal 
slot) is 8 by 1 in., and oriented along the duct length. 

Grid Resolution 

The duct grid is configured to provide the fairest 
comparison among the injectors. Since the injector jet 
shear layers are 3-D and curved in nature, it is not 
possible to construct one duct grid to finely resolve 
these features for each injector case, without an 
extremely large grid. Consequently, the duct portion of 
each grid has near uniform spatial resolution in the 
injector near-field region, and is identical for each 
injector case. It is thought that by employing a nearly 
Cartesian mesh that the development of the turbulent 
flowfield will not be biased towards a given injector 
configuration. 

The injector grids also have nearly consistent spatial 
resolution in all directions, regardless of geometry 
(e.g., each injector face consists of nearly the same 
number of cells. The grid for each injector is “attached” 
to the main duct grid via a mismatched cell face 
boundary. The turbulent boundary layers along the duct 
and injector walls are not resolved in an effort to 
maintain the near-Cartesian look of the grid. This is 
also done because the primary region of mixing is 
thought to be due more to interaction between the 
injection jet and primary stream, as opposed to driven 
by boundary layer effects. 

Using the sequencing feature of the Wind code, a grid 
refinement study was performed (discussed in Results) 
to verify the consistency of the results. The grids could 
be sequenced from coarse to medium to fine with a 


factor of 8 increase in cells (i.e., x2 along I-, J-, and 
K- directions) between each level. Each grid consisted 
of approximately 600,000 points at the finest level. 

Flow Solver 

The flow solver chosen for this study was Wind 
(version 4), which is a merger of three older codes: 
NAST3D [4], NPARC [5], and NXAIR. Wind is a 
product of the NPARC alliance, a partnership between 
the NASA Glenn Research Center and the USAF 
Arnold Engineering Development Center (AEDC) 
dedicated to the establishment of a national, 
applications-oriented flow simulation capability. Wind 
is a node-based, finite-volume, multi-block structured 
code currently administered by AEDC. It permits 
mismatched zonal boundaries, grid sequencing, 
parallel-processing, and multi-species computation, all 
of which are necessary features for the present study. 

A true second-order, spatially accurate (i.e., cell 
stretching is included) numerical scheme is 
implemented to solve the 3-D Reynolds-Averaged 
Navier-Stokes (RANS) equations. The inviscid fluxes 
were computed with Roe flux-difference splitting plus 
limiters to ensure monotonicity in the solution (i.e., to 
satisfy the TVD property). An implicit, 3-D spatial 
factorization scheme, with local time stepping, was 
chosen to drive the solutions to steady state. Mentor’s 
[6] Shear Stress Transport (SST) model was chosen to 
handle the fine-scale turbulence effects; this model is a 
blended combination of k-cb (near walls) and k-e (away 
from walls). 

Numerical Assumptions 

An initial set of runs was conducted using a non- 
reacting 5-species air model. A frozen M = 1.2 
condition was applied at the mainstream duct inflow 
plane, containing pure N 2 . A frozen M = 0.3 condition 
was applied at the injector inflow plane, containing 
pure 0 2 . These multi-species mode runs insured 
constant injector mass flow (recall all injector face 
areas are constant), for the sake of comparison between 
the candidate injectors. Observation and identification 
of the mixing characteristics of the various injector 
configurations was the main goal of this portion of the 
study. Thus, the intent of this portion of the study was 
to determine which configurations generated more 
favorable mixing characteristics, given a constant 
injectant mass flow. An inviscid wall condition was 
applied at the sidewalls of each of the duct grids (to 
enforce the symmetry of a row of injectors). The flow 
conditions at the outflow boundary were extrapolated. 

In a second set of runs for each injector, constant 
stagnation conditions were applied at the injector 
inflow plane. Riemann invariants (one-dimensional) 
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were used to calculate the inflow and outflow 
conditions of the freestream. For these runs, a single- 
species air model was used and effective discharge 
coefficients and induced wall forces were obtained for 
the injector geometries. Table 1 summarizes the applied 
flow conditions. 

RESULTS 

Primary Injector Geometries 

The “primary” set of injector geometries simulated 
includes the circular, diamond, square, and latitudinal 
slot injectors. In each case the massflow through the 
injectors is constant, as explained earlier. The first three 
cases exhibit similar mixing patterns and levels, and are 
distinctly different from the slot result. Figure 3 
displays an injected mass fraction contour plot for the 
circular injector case, showing the flowfield along a 
planar cut through the center of the injector 
(mainstream flow is from left to right). In Fig. 4, a mass 
fraction contour plot is displayed for the same view 
with the latitudinal slot injector. The slot injector fluid 
does not penetrate well, compared to the circular 
injector. Figure 5 displays injected fluid mass fraction 
contour plots for the circular, diamond, square, and 
latitudinal slot injector cases at a plane 50 in. 
downstream of the center of the injector. The counter- 
rotating streamwise vortices carry and mix the injector 
fluid well into the duct flow, causing the increased 
penetration. Interestingly, the diamond injector case 
results in the most complete propulsion of injectant into 
the primary stream, leaving no fluid along the duct 
floor. This more thorough clearing of injected fluid also 
occurs within a shorter distance downstream of the 
injector. Although the diamond injector produces a jet 
that seems to penetrate more quickly, it cannot be 
verified that better mixing is generated. 

A bar chart comparison for the mixing creating by each 
injector is provided in Fig. 6. This “mixing metric” is 
defined herein as the percentage of the primary stream 
flow that contains at least 1 percent injector fluid, along 
a transverse cut 50 in. downstream from injector. Note 
that Fig. 5 illustrates the cuts from which this metric is 
calculated. This metric should be viewed as a 
qualitative indicator of the relative mixing induced by 
each injector. By this measure, the circular, diamond, 
and square injectors produce very similar levels of 
mixing, while the latitudinal slot produces significantly 
less. It is to be noted that latitudinal slot injectors 
are often cited as advantageous for flame-holding 
capability — the current non-reacting simulations cannot 
assess such capabilities. 


The penetration generated by each injector was 
measured using the centroid of injected fluid at various 
downstream positions, following the method proposed 
by Barber et al. [7]. For each constant x plane, fluid 
with an injectant concentration above 1 percent was 
summed to calculate the height of the centroid of 
injected fluid. This centroid penetration is plotted 
against downstream position in Fig. 7. The diamond 
injector consistently displays a slight advantage in 
centroid penetration, followed closely by the circular 
and then square injectors; this trend seems to be 
consistent with the injectant contour plots of the 
primary shapes (see Fig. 5). The latitudinal slot 
produces less than half the centroid penetration of any 
of the other three injector geometries. 

Rectangular Slot Aspect Ratio Study 

The next set of injector geometries simulated involves 
rectangular orifice shapes with varying aspect ratios 
(length/width). Including the latitudinal slot result, there 
are eight rectangular slot cases. Figure 8 displays 
injected fluid mass fraction contour plots for the 
rectangular injector cases at a plane 50 in. downstream 
of the center of the injector. The contour plots show 
improved penetration with increasing aspect ratio. Also, 
the plume appears to spread more, exposing more 
injected fluid to the freestream, allowing a greater 
opportunity for mixing. Both of these assertions are 
supported by the quantitative data. As with the primary 
geometries, the centroid penetration of the injected fluid 
was measured for each injector and plotted in Fig. 9. 
The injectant plume area of each case was calculated 
(using all area with at least 1 percent injectant) and 
plotted in Fig. 10. Clearly, injectant penetration and 
plume spreading increase directly with aspect ratio 
(length downstream). 

The improved penetration is apparently largely due to a 
sustained generation of streamwise vorticity over a 
relatively long streamwise length. And this improved 
vorticity production is likely attributed to a “buffering 
effect.” Specifically, the injector flow closer to the 
trailing edge of a longitudinal rectangular injector is 
better buffered from the axial momentum of the 
freestream by the streamwise vorticities developing 
along the forward portion of the injector. In addition, 
injected flow closer to the leading edge is given some 
upward lift from the trailing edge injector flow blowing 
up from underneath. 

The mixing for the slot injectors is plotted in Fig. 11 
versus aspect ratio. This is the same “mixing metric” 
described earlier. The percentage of mixed fluid also 
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increases with longitudinal aspect ratio. An increase in 
streamwise vorticity production and spreading 
apparently leads to better mixing. As aspect ratio 
increases, a point is reached after which the measured 
mixing appears to approach a limit. A point is also 
reached, in this same range of aspect ratio, after which 
injected fluid is no longer completely cleared from the 
duct floor (refer to Fig. 8). As injector aspect ratio 
increases, the ratio of viscous surface area to fluid 
volume increases (cross-sectional area being held 
constant). This also causes the boundary layer within 
the injector to encompass an increasingly large 
proportion of the injector orifice area. Eventually, the 
viscous effects of increasing this surface area ratio 
would become more detrimental than the benefits of 
increasing the aspect ratio. 

Nonetheless, the longitudinal slot clearly provides 
greater penetration, spreading and mixing than the 
primary geometries. However, for combustor 
application, in which compactness is crucial, a 
longitudinal slot injector may not be practical. 

Effect of Injector Inflow Conditions 

The previous results assumed a fixed, constant 
massflow and flow state for the injector inflow 
boundary. As explained earlier, it was felt that having a 
constant injector-to-mainstream massflow ratio would 
ensure a “fair” comparison of the mixing characteristics 
of each injector. However, this boundary condition is 
not necessarily the most realistic. 

Thus, all cases were recomputed with constant 
stagnation conditions imposed at the injector inflow 
plane. The stagnation pressure and temperature 
correspond to those of the frozen boundary inflow 
cases. This variation allowed for the development of the 
effective discharge coefficient of each case. The 
discharge coefficients are displayed in Table 2. The 
circular, diamond, and square injector cases all have 
very similar discharge coefficients and there is no 
indication of any one trending toward particularly 
greater or lesser values. Nonetheless, all three of these 
injectors produce a considerably larger discharge 
coefficient than the latitudinal slot injector. As for the 
rectangular slot injectors, a clear trend is present. There 
is a direct relationship between discharge coefficient 
and aspect ratio, though the discharge appears to 
approach a maximum at the case with largest aspect 
ratio (8:1). This maximum discharge phenomenon is, 
no doubt, caused by the increase in viscous effects 
associated with increasing the ratio of injector surface 
area to fluid volume, which occurs as a rectangle’s 
aspect ratio increases (as discussed earlier). 


Induced Wall Force Study 

The data from this study has been analyzed with a 
primary focus on mixing applications, since mixing 
drives combustion and perhaps noise reduction. 
However, thrust-vectoring of low speed supersonic 
flows is an application of great interest. An integrated 
non-dimensional vertical wall force was computed for 
each injector case (shown in Table 3). This wall force 
was computed for the 20-inch long section of the duct 
floor that contains the injector orifice. 

The whole duct floor would have been used in the 
computation, but a difficulty with one of the boundary 
conditions arose in the solutions. In Fig. 12, a Mach 
contour plot is shown of the diamond injector case. An 
expansion wave is generated on the back lip of the 
injector orifice. This expansion wave propagates itself 
through the solution and then reflects off of the 
boundary at the duct ceiling. When the wave reflects off 
of the ceiling boundary condition, it propagates itself 
back into the downstream region of the solution and 
corrupts the pressure data computed on the floor of the 
duct downstream of x = 40 in. Consequently, only the 
floor region between x = 20 and x = 40 in. was used for 
wall force calculation. 

The circle, diamond, and square injectors generate very 
similar wall forces, considerably greater than that of the 
latitudinal slot. As for the rectangular injectors, the 
induced wall force seems to increase with aspect ratio 
and approach something of a limit near the largest 
aspect ratio of 8:1. It is possible that this phenomenon 
is simply related to the increase in injector mass flow 
from the longer aspect ratio rectangular slot. 

Grid Convergence Study 

For all the cases in this study, grid independent results 
could not be obtained without an excessively large grid. 
Consequently, the data acquired is primarily of 
qualitative value. As described earlier, coarse, medium, 
and fine grid solutions were obtained for each case. 
Whereas accurate quantities could not be computed, 
comparisons were made based on all three grid levels to 
ensure the constancy of the observed trends. The trends 
observed in this study — mixing behavior, induced wall 
forces, and discharge coefficients — are all consistent at 
all grid resolutions. Both medium and fine grid 
calculations of discharge coefficient, lift coefficient, and 
mixing, for all cases, can be found in Tables 2, 3, and 4, 
respectively. For future study, it may be possible to apply 
an adaptive grid strategy to obtain grid independent 
results without excessive computational cost. 
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SUMMARY 

The circular, diamond, and square injectors provide 
deeper injected fluid penetration than does the wide slot 
injector. This appears to be due mainly to the lifting 
effects of the increased streamwise vorticity these 
geometries generate. Fluid injected through the 
diamond injector penetrates more rapidly than that 
through the circular and square injectors and is more 
thoroughly cleared from the duct floor; this is consistent 
with the findings of others for transverse injection 
through diamond injectors at similar pressure ratios. 
Though the diamond injector generates more rapid 
penetration, it cannot be determined, from this study, 
whether or not it produces better mixing. 

For rectangular injectors, the penetration depth and 
measured mixing increase with aspect ratio (L/W). 
Again, this is largely due to a sustained generation of 
streamwise vorticity over a relatively long streamwise 
length. This improved vorticity production is likely 
attributed to a buffering effect, in which the aft portion 
of the injected flow is shielded from the mainstream by 
the streamwise vortices, which originate along the 
forward portion of the injector. The measured discharge 
coefficient and induced wall force also increase with 
aspect ratio, though the increase in measured wall force 
may simply be a function of increased discharge. 
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Table 1: Summary of freestream and in 

ector flow conditions 


Freestream 
(all runs) 

Fixed static injection 
(1st set) 

Fixed stagnation injection 
(2nd set) 

M 

1.2 

0.3 

— 

P (psi) 

14.7 

29.4 

31.3 

T(°R) 

530.0 

670.6 

682.6 


Table 2: Discharge coefficients 


Configuration 

Medium grid 

Fine grid 

Circle 

0.460 

0.528 

Diamond 

0.488 

0.524 

Latitudinal slot 

0.335 

0.287 

Rectangle 1:3 

0.333 

0.333 

Rectangle 1:2 

0.418 

0.443 

Rectangle 1:1 (square) 

0.464 

0.501 

Rectangle 2:1 

0.517 

0.553 

Rectangle 3:1 

0.568 

0.636 

Rectangle 5:1 

0.608 

0.712 

Rectangle 8:1 

0.605 

0.719 
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Table 3: Lift coefficients 


Configuration 

Medium grid 

Fine grid 

No injection 

0 

0 

Circle 

-0.1194 

-0.1197 

Diamond 

-0.1213 

-0.1251 

Latitudinal slot 

-0.09492 

-0.07090 

Rectangle 1:3 

-0.07529 

-0.07364 

Rectangle 1:2 

-0.09975 

-0.1025 

Rectangle 1:1 (square) 

-0.1244 

-0.1179 

Rectangle 2:1 

-0.1377 

-0.1344 

Rectangle 3:1 

-0.1479 

-0.1563 

Rectangle 5:1 

-0.1660 

-0.1692 

Rectangle 8:1 

-0.1527 

-0.1744 


Table 4: Mixed flow percentage 


Configuration 

Medium grid 

Fine grid 

Circle 

17.670 

13.378 

Diamond 

17.032 

13.624 

Latitudinal slot 

11.814 

8.822 

Rectangle 1:3 

14.946 

12.184 

Rectangle 1:2 

15.873 

13.346 

Rectangle 1:1 (square) 

17.193 

14.341 

Rectangle 2:1 

19.177 

15.285 

Rectangle 3 : 1 

20.260 

16.246 

Rectangle 5 : 1 

22.439 

18.649 

Rectangle 8 : 1 

27.839 

20.539 



Freestream 
M = 1.2 


Injector 
M k 0.3 


Figure 1: Supersonic duet with transverse injector 
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Figure 2: Injector orifice geometries 
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Figure 4: Mass fraction contour plot for duct with 1:8 slot injector 
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Figure 5: Mass fraction contour plots for injector ducts 50 in. downstream of injector 
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Figure 6: Mixed flow percentage at 1 percent minimum injectant mass fraction tolerance 
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Figure 7: Injectant centroid penetration 
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L/W = 0.125 


L/W = 0.36 


L/W = 0.59 


L/W = 1 



L/W = 1.7 


L/W = 2.8 


L/W = 4.8 
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Figure 8: Injectant mass fraction contours of rectangular injector ducts 50 in. downstream of injector 



Figure 9: Centroid penetration (rectangular injectors) 
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Figure 10: Jet area (rectangular injectors) 
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Figure 11: Percent flow mixed versus aspect ratio for rectangular injectors 
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Figure 12: Mach contour plot 
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